Robustness of quantum correlations against decoherence 
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We study dynamics of nonclassical correlations by exactly solving a model consisting of two atomic 
qubits with spontaneous emission. We find that the nonclassical correlations defined by different 
measures give different qualitative characterizations of those correlations. The relative behaviors 
of those correlation measures are presented explicitly for various quantum states in the two-qubit 
atomic system. In particular, we find that the robustness of quantum correlations can be greatly 
enhanced by performing appropriate local unitary operations on the initial state of the system. 
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I. INTRODUCTION 

Entanglement plays a key role and is a resource 
in quantum information processing (QIP). In the last 
decades, much efforts have been put into study it in vari- 
ous physical systems, see for example Refs. [I-Q- While 
being widely considered to be something unique without 
classical counterpart, it is recently realized that entangle- 
ment reflects only one aspect of quantum correlation of 
a quantum state, there exist some other nonclassical cor- 
relations (6l-[l0|. One example is a model of QIP in Ref. 
[ll|, which outperforms its classical counterpart while 
without any entanglement. It is proposed that quantum 
discord (QD) ^6i|, a different measure of quantum corre- 
lation, might be responsible for the power of this model 
p^ . Thus, entanglement and QD, and additionally some 
other measures all are important in describing nonclassi- 
cal correlations for a quantum system. 

It is well acceptable that decoherence will cause de- 
creasing of entanglement which may induce failure of the 
algorithms and various protocols of QIP. On the other 
hand, it is also of fundamental interest to study the 
behavior of nonclassical correlations with decoherence. 
Recently, there are some studies in evaluating the rela- 
tionships between nonclassical correlations quantified by 
QD and entanglement quantified by concurrence or en- 
tanglement of formation (EoF) under the action of de- 
coherence [T3l - [2^ . Also the role of QD in identifying 
critical points of quantum phase transitions are studied 
[25l - l28j . For practical use, one may hope that the non- 
classical correlations which are crucial for QIP can be 
maintained for sufficiently long time so that the designed 
tasks can be fulfilled. However, the unavoidable inter- 
action of a realistic system with its surroundings always 
makes it decaying with time. Particularly, under certain 
circumstances, the entanglement of a bipartite state can 
even terminate abruptly in a finite time, a phenomenon 
termed as entanglement sudden death (ESD) by Yu and 
Eberly For nonclaccical correlations as measured by 
QD and other similar quantities, however, the sudden 
death phenomenon has not been observed (this may be 



good for QIP independent of entanglement), but their 
values may decay with time and become zero at discrete 
instants of time Thus it is of great importance to 

find ways to preserve them when the system is in contact 
with an environment which leads to decoherence. 

In this paper, we offer a comparative study of the rela- 
tionships between entanglement and different measures 
of nonclassical correlations for a certain family of two- 
qubit states. The system we considered consists of two 
atomic qubits with spontaneous emission [2^. We will 
show that even in cases where entanglement disappears, 
nonclassical correlations persist in the whole time region. 
In this sense, they are more robust against decoherence 
than entanglement, implying that quantum algorithms 
based only on them may be more robust than those based 
on entanglement. Moreover, we will show by explicit ex- 
amples that robustness of quantum correlations can be 
greatly enhanced by performing local unitary operations 
to the initial states, though it is well-known that those 
measures of correlations themselves are invariant under 
local unitary operations 

This paper is organized as follows. In Section II we 
recall some measures of correlations for a bipartite state. 
We will adopt the EoF as the measure of entanglement, 
and other nonclassical correlations will be quantified by 
the measurement induced disturbance (MID), QD and 
the geometric measure of QD (GMQD). In Section HI, we 
introduce the model and compare characteristics of var- 
ious correlation measures for a class of two-qubit states. 
We will also evaluate robustness of quantum correlations 
for different local unitary equivalent states. Finally, Sec- 
tion IV is devoted to a summary. 



II. MEASURES OF QUANTUM 
CORRELATIONS 



The basic and crucial problem for studying dynamics of 
correlations of a composite system is to characterize and 
quantify correlation from different point of views in that 
system. The motivation of this study is to find whether 
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those physical systems can be a platform for various tasks 
of QIP. In this section, we will first list several measures of 
quantum correlations for the bipartite states, and those 
measures will be applied in our study of the two-qubit 
atomic system. 

First, we still consider entanglement as one important 
quantum correlation. In order to study the entangle- 
ment dynamics for any two-qubit states, we adopt EoF [1] 
which is a well-accepted measure of entanglement. The 
definition is as follows, 

where H{t) = — Tlog2 r — (1 — r) log2(l — r) is the bi- 
nary Shannon entropy, and C is the time-dependent con- 
currence 0, which takes the form C = max{0,Ai — 
A2 — A3 — A4}, where Ai (« = 1,2,3,4) are the square 
roots of the eigenvalues of the spin-flipped operator R — 
p{a'^ ® a^)p*{a^ (g) cr^) arranged in decreasing order, with 
p* being the complex conjugation of p in the standard 
basis, and is one Pauli matrix. 

Besides entanglement, let us consider some other quan- 
tum correlations. The so-called MID, proposed recently 
by Luo in Ref. 0, is defined as the difference of two 
quantum mutual information respectively of a given state 
p shared by two parties a and b and the corresponding 
post-measurement state Il{p), 

MID(p) = /(p)-/(n(p)), (2) 

where the mutual information, I{p) — S{p"-) + S{p'') — 
S{p), measures the total correlation, including both clas- 
sical and quantum, for a bipartite state p. Here S{p) = 
— tr(/5log2 yo) denotes the von Neumann entropy, with 
and p^ being the reduced density matrix of p by tracing 
out b and a, respectively. While quantum mutual infor- 
mation of the state n(p), /(n(p)), quantifies the classical 
correlation in p, with U{p) = (Hf (g) n^^)/3(nf (g) Hj), 
where the measurement is induced by the spectral res- 
olutions of the reduced states = Eipfll" and p^ = 
Tijp^H'j. In fact, Il{p) is a classical state (a state is re- 
ferred to as classical if it is invariant under local von 
Neumann measurement) for any complete set of projec- 
tive measurements {Hf} and {IlJ}- Il(yo) defined above is 
actually the closest classical state to p since this kind of 
measurement leaves the marginal information invariant 
and is in a certain sense the least disturbing |7| . 

Another well- accepted measure of nonclassical corre- 
lation is the QD 6]. It was defined by the discrepancy 
between quantum mutual information and the classical 
aspect of correlation which was defined as the maxi- 
mum information about one subsystem that can be ob- 
tained by performing a measurement on the other sub- 
system. If we restrict ourselves to the projective mea- 
surements performed locally on subsystem a described 
by a complete set of orthogonal projectors {11^}, then 
the quantum state after measurement changes to pb|fc = 
(Ilfc (g) T)p{Ilk <E) I) /pk, where I is the identity operator for 



subsystem b, and pk = tr[(nfc (E) I)p{Ilk g) I)] is the prob- 
ability for obtaining the measurement outcome k on a. 
The classical correlation can then be obtained by max- 
imizing Jip\{Uk}) = S{p'') - Sip\{Uk}) over all {n^}, 
where 5(/9|{nfe}) = T,kPkS{pb\k) is a generalization of 
the classical conditional entropy of subsystem b. Explic- 
itly, the QD is defined as the minimal difference between 
I{p) and J(p|{nfc}) as 

QD(p)=/(p)- sup J(p|{nfc}), (3) 
{n^} 

where the supremum is taken over the complete set of 
{Ilfc}. The intuitive meaning of QD thus may be inter- 
preted as the minimal loss of correlations due to measure- 
ment. It vanishes for states with only classical correlation 
and survives for states with quantum correlation. 

Since analytical expressions of QD are achievable only 
for certain special classes of states [2l|, [22j , Dakic et al. 
[l] proposed the GMQD. They use the square of the 
Hilbert-Schmidt norm as the distance between two quan- 
tum states, and the nearest distance between p and all 
of the zero-discord states is interpreted as a measure 
of quantum correlation, which is defined as 



GMQD(p)- min ||p-xf, (4) 

where the geometric quantity || p — x IP= tr(p — x)"^ ■ 

By noting that any two-quit state p can be represented 
asp = + X ■ a ^I + I(Ey- a + J2tj=i '^^i'^i '^i)' 

Dakic et al. derived an explicit formula of GMQD, which 
is given by ^ 

GMQD(p) = i (II X ||2 + II i? ||2 , (5) 

where the scalar product a-a = "^jf j with a — x,y. 

Moreover, ||x|p = J2^=ixh II-^IP = tr(i?^i?), and /cmax 
is the largest eigenvalue of the matrix K = xx^ + RBF , 
where the superscript T denotes transpose of vectors or 
matrices. Since this measure gives analytic results with- 
out restriction on the form of p, it turns out to be a con- 
venient tool for analyzing quantum correlation dynamics 
from a geometric perspective. 

Correlation measures (i.e., EoF, MID and QD) listed 
above assume equal values for bipartite pure states, but 
this is not the case for mixed states. Particularly, the QD 
and GMQD are not symmetric quantities with respect to 
the measurements performed on subsystem a or 6, and 
the state is said to be completely classically correlated 
only when both of them reach zero. In the subsequent 
discussions, however, we will restrict ourselves to the sit- 
uation that the states are symmetric under exchange of 
subsystems. 
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III. CORRELATION DYNAMICS UNDER 
DECOHERENCE 

We are interested in revealing difference between vari- 
ous correlation measures under the action of noisy envi- 
ronments. For this purpose, we consider a system con- 
sists of two identical atoms with spontaneous emission, 
and having lower and upper levels \gi) and \ei) {i = 1, 2) 
separated by the energy gap huj, with cj being the tran- 
sition frequency. Also, we assume this system is coupled 
to a multimode vacuum electromagnetic field, under the 
influence of which its time evolution is governed by the 
following master equation [2^ [s^] 



dp 
dt 
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1=1 
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1 2 

~-J2j,,{2S7pS+-{S+S;,p}), 



(6) 



where Sf^ are the raising and lowering operators while S'f 
is the energy operator, jij = 7 (« = j) are the sponta- 
neous emission rates of the atoms caused by their direct 
interaction with the vacuum field. Moreover, 7^^ and fiy 
(i j) describe the collective damping and the dipole- 
dipole interaction potential, respectively. They both de- 
pend on the interatomic distance Vij = \rj — r^j and take 
the form 



sin(fcrij) ^ cos(krij) 

kvij 
sin(fcr,ij) 



sin(fcrij) 



cos{krij) cos{krij) 
[kr^jY kr.,j 



(7) 



where we have assumed that the atomic dipole moments 
for the two atoms take the same value (i.e., fJ.i = IJ.2 = p) 
and are polarized in the direction perpendicular to the 
interatomic axis. And, k = 2-k /\ is the wave vector with 
A being the atomic resonant wavelength. 

We first consider the Bell-like initial state as follows 



I*) = a|eie2) + Vl-o^l. 91.92), 



(8) 



for which the nonzero elements of p{t) in the standard 
basis {\e1e2), \e1g2), \g1e2), \g1g2)} are of the form 



Pii{t) 
Pii{t) 

P22,33{t) 
P23,32{t) 



Pliit) 



71+2* 



e - '-J + 02 [e '^12* - 



-27*1 



-27*1 



a2 e 



-7l2* 



-27*1 



(9) 



where 7^^ = 7 ± 712, ai.2 — c?'l\2l'^l\2- one can see, 
the X structure of p{€) is maintained during the evolution, 
which greatly facilitates the following analysis. 

Fig. 1 shows the contour plots of the quantum correla- 
tions versus 7^ and o? in the short time region with the 
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FIG. 1: (Color online) Contour plots of different correlation 
measures versus and cP' for the initial state where the 
interatomic distance for these plots is 7-12 = 0.6737A. 



interatomic distance r\2 = 0.6737A, at which 712 reaches 
its minimal value. Clearly, while the ESD happens after 
a finite time, the MID, QD and GMQD maintain during 
the whole time region. This points to a fact that the non- 
classical correlations may be more resistant to external 
perturbations than that of entanglement. Moreover, the 
GMQD shows a revival for large values of a^, but dur- 
ing the same time region the other correlation measures 
are always decreased. Thus similar to the relativity of 
different entanglement measures (see [3]| and references 
therein), the nonclassical correlation measures may also 
impose different orderings of quantum states. 

For quantum states of Eq. (9), the concurrence can be 
obtained analytically as C{p) = 2 max{0, Ci, C2}, where 
Ci = |pi4| - VP22P33 and C2 = \p23\ - VPii/°44- In 
the short time region, |pi4p > P22P33, thus the concur- 
rence as well as EoF are continuous functions of ^t. In 
the long time region, I/O23P > Pii/344, and the entangle- 
ment may experiences a very weak revival after a finite 
time interval of its complete disappearance. Similarly, for 
p{t) of Eq. (9) we have Mm{p) = S(n{p)) ~ S{p), with 
n(p) = diag{pii, P22, P33, /344}, thus it is understandable 
that the MID also behaves as continuous functions of -ft 
and a^. For the QD and GMQD, as denoted by the blue 
and green curves in Fig. 2, they present sudden changes 
at discrete instants of time or which are evidenced by 
the presence of kinks. In general, these sudden change 
behaviors are caused by the optimal procedure for choos- 
ing the measurement operators, and we will discuss this 
issue in detail in the following. 

From Fig. 2 one can also observe that the MID, QD 
and GMQD behave as non-monotonic decreasing func- 
tions of jt and . Their values may be increased or de- 
creased during certain parameter regions. Particularly, it 
is worthwhile to note that there are states having differ- 
ent orderings induced by different correlation measures. 
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FIG. 2: (Color online) Quantum correlations versus 'yt with 
= 0.9 (top) and versus o? with 'yt = 0.35 (bottom) for 
the initial state |^) with ri2 = 0.6737A. The insets show the 
cases when 7i £ [1, 6] and £ [0.97, 1], where the curves are 
cut to better visualize of the revival phenomenon, and the 
vertical lines indicate the critical 'yt at which the correlations 
turn out to be increased or decreased. 

For example, for the parameters chosen in Fig. 2 there 
exists states pi and p2 such that MID(pi) > MID(p2), 
QD(pi) > QD(p2), while GMQD(pi) < GMQD(p2) 
when 7^ e [0.639,0.888], i.e., the decrease of MID and 
QD are accompanied by the increase of GMQD. More- 
over, from the insets of Fig. 2, other states differently 
ordered by the nonclassical correlation measures are read- 
ily distinguished, including those during the scaled time 
regions -ft e [2.304,2.604] and -ft G [3.672,4.164], for 
which we have MID(pi) > MlD{p2) and GMQD(pi) > 
GMQD(p2), while QD(pi) < QD(p2)- Thus we see that 
there are no simple dominance relations between MID, 
QD and GMQD for general cases. They are different not 
only quantitatively but also qualitatively. 

Now we evaluate analytically the sudden change be- 
haviors of the QD. For p{t) in Eq. (9), the infimum of 
the conditional entropy can be determined as 

inf 5(p|{na) = i/(r), (10) 

where the optimal angle 9 related to the projective mea- 
surement \s 9 = {m + 0.5)7r with m G Z, and H(t) is the 
Shannon entropy functional, with 

1 - v/[l - 2(pn + P33)]^ + 4(|pi4|^ + \P2Z? + 5) 

^= 2 ' 

(11) 

with (5 = 2a(l-a2)i/2e-'^V23COs(2(^-2wi). Clearly, the 
infimum of 5(/3|{nfc}) is obtained whenever (/) = ki: + ujt 



FIG. 3: (Color online) QD versus yt with ri2 = 0.6737A, 
= 0.9 (top) and versus ri2/X with •yt = 0.35, = 0.8 
(bottom) for the initial state l^*). 

if ap23 > 0, and (j) = {k + 0.5)tt + ujt (fee Z) if ap23 < 0. 
Thus we see that the sudden change behaviors of QD here 
is caused by the optimization procedure for choosing the 
optimal measurement angle (j) or equivalently, the optimal 
projective measurement operators {11^} because for the 
two-qubit case 11^ can be expressed aslli = {I + n- a)/2 
and 112 = I — Hi , with n = (sin 6 cos (j), sin 6 sin 0, cos 6)'^ 
being a unit vector in R'^. Here it should be emphasized 
that the angle that minimizes the QD is time depen- 
dent, which is different from the previous results with 
real elements of the density matrix (Tsl - flTj . As is shown 
evidently in the top panel of Fig. 3 with fixed ri2 and 
a^, the optimal measurement angle (f> assumes the value 
(j) — (k + 0.5)7r -f- ujt in the short time region, while in the 
long time region it is (/> = kir + Lut. For fixed jt and a^, as 
exposed in the bottom panel of Fig. 3, the two different 
choices of 4> yields the QD alternatively. 

When considering the GMQD, from Eq. (5) it is read- 
ily to see that it can also be expressed as GMQD(/9) = 
(tiK - fcinax)/4, where tri^ = T,ik, and /cmax = 
maxj/ci, fc2, fcs}. For quantum states with the X struc- 
ture (of course, including those of the form of Eq. (9)) 
one can obtain the eigenvalues of K analytically as 

ki^2 = 4(|pi4|±|p23|)', 
4 

h = 2^p2^-4(piiP33+p22P44)- (12) 

The ki are in arbitrary order and cannot be ordered 
by magnitude unless the parameter values are known, 
thus /cmax may be changed with the variation of the pa- 
rameters involved, which causes the discontinuities of the 
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FIG. 4: (Color online) GMQD versus with = 0.9 (top) 
and versus o? with -yt = 0.35 (bottom) for the initial state 
with the interatomic distance r\2 = 0.6737A. 



GMQD. Two exemplified plots for certain specified sys- 
tem parameters are presented in Fig. 4, from which one 
can see that for fixed — 0.9, the GMQD is obtained 
when fcmax = ^3 in the short and long time regions, while 
in the middle time region, it is obtained when fcmax = fci ■ 
For fixed = 0.35, we have fcmax = ^3 for small or large 
values of a^, and fcmax = for middle values of . 

In Ref. Dakic et al. showed that /cmax can be writ- 
ten as fcmax — max|g>=i| Ke, and the sudden change in 
decay rates of GMQD corresponds to the sudden change 
of the optimized e — (sin cos </>, sin 9 sin 0, cos 9)^ in R'^. 
For the two-qubit case, the closest zero-discord state 
X " Yij,^iPk^k ® Pk is obtained when Hi = (I -|- e • (t)/2 
and 112 = I — Hi, with e being the eigenvector of K with 
the largest eigenvalue. The eigenvectors correspond to 
the three eigenvalues of Eq. (12) are 



I , \ A Re P14P23 T P14P23 „ 

V Im(pi4P23j 

IV3) = (0,0,1)^, 



(13) 



where ai 2 are normalization constants, e assumes 
or l^/'a) with the variation of the involved system param- 
eters. If e = 1-01) then we have cos6' = or equivalently, 
9 = (n + 0.5)7r with n e Z. If e = \'4'z) , however, we have 
cos 9 = ±1 or = riTT with n ^ 1. The optimal cf) can 
also be derived exactly, and here we do not care about it 
for the sudden change behaviors of GMQD have already 
been reflected by the discontinuity of the optimal 9. 

Notice that although both the optimal unit vectors n 
for the QD and e for the GMQD are related to the pro- 
jective operators, they are not the same for general case, 
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FIG. 5: (Color online) Quantum correlations versus ri2/X for 
the initial state |^) with — 0.8 and jt = 0.35. 




FIG. 6: (Color onhne) Quantum correlations versus -yt for 
the initial state p*(0) = p|vE'){*| + (1 -p)I/4 with = 1/2, 
p = 0.65 and ri2 = 0.6737A (red curves). The cyan regions 
show the extent to which the correlations can be adjusted by 
performing operations {I x V} to the second qubit. 



which is reflected in the disagreement of their individual 
sudden changes in the decay rates. 

In Fig. 5 we show the influence of rxijX on correla- 
tion dynamics for the initial state \^) with o? — 0.8 and 
— 0.35. It is seen that all the correlation measures be- 
have as damped oscillations with increasing ri2/A, and 
the dependence of QD and GMQD on ri2/A are more sen- 
sitive than those of EoF and MID. Also one can note that 
there exists states p\ and p2 such that QD(pi) > QD(/32), 
GMQD(pi) > GMQD(p2) and MID(pi) < MID(p2)- 
Thus we see again that the states have different order- 
ings induced by different measures of nonclassical corre- 
lations. 

Finally, we show that the robustness of quantum corre- 
lations may be enhanced by performing local unitary op- 
erations to the initial state p(0). As we know, the corre- 
lation measures (EoF, MID, QD and GMQD) are locally 
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unitary invariant; that is, p{0) and {U (^V)p{0){W (^V^) 
have completely the same values of correlations for any 
unitary operators U and V acting on parties a and b. But 
this does not mean that their robustness against decoher- 
ence are the same. In fact, for any initial state p(0), there 
exist an optimal local unitary transformation [/opt 'Si Vopt 
which can enhance the robustness of the correlations to a 
certain greatest extent. However, in general, the optimal 
Uopt S Vopt is determined by the decoherence mechanism 
as well as the explicit form of p{0), and it is difficult to 
obtain an analytical expression of it. 

For simplicity, we consider the initial state as a Werner- 
like state p*(0) = + (1 -p)I/4. And we restrict 
ourselves to the complete set of locally unitary transfor- 
mations {I X V} (a subset of {U €5 V}). Density matrix 
for p*(0) can be obtained analytically [30], and the 
dependence of EoF, MID, QD as well as GMQD on the el- 
ements of p* (t) are the same as those for the initial state 
l^*). But the parameter appeared in Eq. (11) changes 
to (5 = 2|/9i4||p23| cos(2(/) — wi — W2), where cosa;i_2 = 
Re(pi4,23)/|pi4,23| and sinwi,2 = -Im(pi4,23)/|pi4,23|- 
For general p{0) S {(I x V)p^{0){I x V^)}, however, the 
density matrix p{t) as well as the correlation measures 
can only be obtained numerically since the X structure 
of it may be destroyed under the transformation {I x V}. 
By expressing V in terms of the Pauli operators as 
V — sma{(T^ cos 7 -I- cr^ sin 7) -I- cos a{a^ cos /? — i(T° sin /3) 
with a G [0, 7r/2] and l3,"f [0, 2tt], we perform numerical 
calculations with p — 0.65, — 1/2, ri2 — 0.6737A and 
show the results in Fig. 6. It is shown that the extent of 
different correlation measures at finite time t > can be 
adjusted among values bounded by the cyan regions by 
performing operations {I x V} to the initial state. Par- 
ticularly, the sudden death phenomenon of entanglement 
may be avoided. Thus we see that different locally uni- 
tary equivalent states show different robustness against 
decoherence induced by the spontaneous emission, and 
this provides a way for long-time preservation of quan- 
tum correlations in this system. 

IV. CONCLUSIONS 

In conclusion, we evaluated systematically the quan- 
tum correlations quantified by the EoF, MID, QD and 



GMQD by exactly solving a model with the two qubits 
interacting with a multimode radiation field. We showed 
that the dynamics for these quantities may be rather dif- 
ferent. Particularly, while entanglement experiences sud- 
den death, MID, QD and GMQD persist for an infinite 
time, i.e., they are more robust than that of entanglement 
against decoherence, which is important for QIP tasks 
which deos not based on entanglement. Moreover, the 
dissipative process of spontaneous emission may lead to 
a generation of states manifesting the relativity of differ- 
ent correlation measures. In general, different measures 
of nonclassical correlations are incomparable since their 
behaviors may be qualitatively different. 

We also demonstrated the disagreement between QD 
and GMQD on reflecting the sudden change behaviors of 
nonclassical correlations of a given quantum state. The 
origin of this phenomenon is the different definitions of 
discord based on different measures of distance between 
two states. While the QD is based on the von Neumann 
entropy, the GMQD is based on the Hilbert-Schmidt 
norm, and it is this difference that yields different op- 
timal measurement operators. 

Since it is generally accepted that quantum correla- 
tions is crucial for QIP, it is important to maintain them 
as long as possible. Here we showed that although local 
unitary operations cannot change extents of quantum 
correlations of a state, it can greatly enhance the robust- 
ness of a system against decoherence. This provides a 
possible way for long-time preservation of the quantum 
correlations by local unitary operations. 
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